Computing the Casimir energy using the point-matching method 
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We use a point-matching approach to numerically compute the Casimir interaction energy for 
a two perfect-conductor waveguide of arbitrary section. We present the method and describe the 
procedure used to obtain the numerical results. At first, our technique is tested for geometries with 
known solutions, such as concentric and eccentric cylinders. Then, we apply the point-matching 
technique to compute the Casimir interaction energy for new geometries such as concentric corru- 
gated cylinders and cylinders inside conductors with focal lines. 
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I. INTRODUCTION 



One of the most clear macroscopic manifestations of 
quantum mechanics is the Casimir force a tiny force, 
relevant at short distances, which appears on uncharged 
bodies due to changes in the zero-point energy associ- 
ated to quantum vacuum fluctuations. Quantum effects 
hkc Casimir forces have become increasingly important 
as electronic and mechanical systems on the nanome- 
ter scale become more prevalent. In the last 10 years, 
Casimir-force measurements have been made in a wide 
variety of experiments. For example, plate-plate 0, 
sphere-plate [3[ and crossed cylinders [J] configurations. 
There are also proposals for measuring the Casimir force 
in the cylinder-plane configuration [5|. A recent exper- 
iment ^Q] implies that Casimir interactions arc relevant 
in the fabrication of commercially MEMs devices. Tech- 
niques for predicting Casimir forces in general geome- 
tries are clearly needed if theory is to keep pace with the 
flourishing of experimental Casimir data and the design 
challenges of future MEMs devices. 

In this framework, accurate modeling of Casimir forces 
in general geometries and for arbitrary electromagnetic 
properties of the uncharged bodies becomes very impor- 
tant. Shape and geometry can strongly influence the 
Casimir interactions. Several theoretical techniques have 
been developed in order to understand the geometric de- 
pendence of the Casimir force. These include the use of 
the argument theorem to perform explicitly the sum over 
modes M, H, H, , semiclassical and optical approxima- 
tions methods based on functional integrals and 
scattering theory [l^. Many of these approaches have 
a common root in the multiple scattering theory devel- 
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oped in the seventies [IJI (see also [1^1 for an updated 
review and applications to semitransparent bodies), and 
the evolution in the computational power allowed a pre- 
cise numerical evaluation that involves, in general, the 
computation of determinants of infinite matrices. 

There are also full numerical approaches, as the world- 
line numerics [TB], that has been applied to scalar fields 
satisfying Dirichlet boundary conditions, or finite differ- 
ence methods that evaluate the Casimir energy from the 
two-point function of the electromagnetic field As 
a consequence of this theoretical activity, we now have 
exact results for a variety of geometries that involve per- 
fectly conducting shells: cylinder and sphere in front 
of a plane ITsl. Il9l. [20t , eccentric cylinders [13, HH , 
two-spheres jl9l |22| , surfaces with periodic corrugations 
[23! , Casimir pistons [l^l , multibody interactions involv- 
ing plates and cylinders [H] or sauar es [l7| . objects of 
spheroidal or nearly-spheroidal shape [26| . etc. Some of 
these methods also apply to the case of imperfect mirrors. 

The Casimir force between two conductor bodies with 
simple shapes is in general attractive and monotonically 
decreasing with the separation among conductors. Thus, 
one might ask whether complex geometries might give 
rise to unexpected phenomena, such as non-monotonic 
forces. For more complicated geometries, however, full 
numerical calculations become extremely difficult and 
therefore it is worth to analyze alternative approaches. 

In this paper, we explore a new numerical approach to 
compute the Casimir interaction energy in a two conduc- 
tor waveguide for arbitrary cross section. The idea is to 
combine a well known method for computing eigenvalues 
of the Helmholtz equation, the point-matching technique 
[131, with the argument theorem to compute explicitly 
the sum over eigenfrequcncies. Our technique is first 
tested for geometries with known solutions and then is 
applied to new geometries. Thus, in Section |ll] wc will 
describe the new approach. In Section lllll we validate our 
method against numerical results [lol . l2l| obtained by us 
in previous works. In the rest of the paper, we use our 
numerical approach to compute the Casimir interaction 
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energy for more general cross sections of the waveguide, 
as corrugated cylinders in Section |TVl In Section IV Al 
and IV B( we explore the behaviour of the Casimir inter- 
action energy between outer conductors that have focal 
lines and a cylindrical inner conductor. Concretely, Sec- 
tion IV Al shows the Casimir interaction energy when a 
cylinder with circular section is placed inside an outer 
cylinder with elliptical section. In Section IV B I the outer 
conductor is a waveguide with parabolic section. Finally, 
in Section IVTl we expose our final remarks. 



II. POINT-MATCHING NUMERICAL 
APPROACH 

A two-conductor waveguide presents an interesting 
setup for the application of the point-matching technique. 
This technique has been widely used to solve eigenvalue 
problems in many areas of engineering science [27|. The 
boundary conditions are imposed at a finite number of 
points around the periphery of both conductors. Un- 
der this assumption, and for a solution proportional to 
g-iwt^ the partial differential equation of the problem can 
be reduced to a system of linear algebraic equations. The 
determinant associated to this system vanishes for some 
values of w, the eigenfrequencies of the system, that can 
be determined in this way by searching numerically the 
zeros. In order to obtain the Casimir energy, instead of 
computing each eigenvalue, it is more convenient to use 
the argument theorem to perform the sum over all eigen- 
values. 




FIG. 1: A two-conductor waveguide in which one conductor 
encloses the other and each has arbitrary cross section. 

Let us consider a general geometry with translational 
invariance along the z— axis (as for example very long 
and parallel conducting shells of arbitrary sections). We 
will bear in mind the situation in which one conductor en- 
closes the other, as shown in Fig. [Tl although the method 
could be applied to more general cases. The Casimir in- 
teraction energy for a system like this, composed of two 
conducting shells, can be written as 

P 



where Wp are the eigenfrequencies of the electromag- 
netic field satisfying perfect conductor boundary con- 
ditions on the surfaces of the conductors, and Wp are 
those corresponding to a situation in which the exter- 
nal conducting shell is very large. Throughout this pa- 
per we use units h = c = I. The subindex p denotes 
the set of quantum numbers associated to each eigenfre- 
quency. Introducing a cutoff for high frequency modes 
£^12(0-) = 5 J2pi'^~'^"''"'^p ~ e~'^'^'"Wp), the Casimir inter- 
action energy £'12 is the limit of i?i2(o') as ct ^ 0. For 
simplicity we choose an exponential cutoff, although the 
explicit form is not relevant. 

The transverse electric (TE) and transverse magnetic 
(TM) modes can be described in terms of two scalar 
fields, with adequate boundary conditions. In cylindrical 
coordinates, the modes of each scalar field will be of the 
form hn,kAt^r,9,z) ^ e(-™"''==*+*'=-^)i?„(r, 61), where the 
eigenfrequencies are Wn.k, = + A^, and A„ are the 
eigenvalues of the two dimensional Laplacian 

The set of quantum numbers p is given by (n, fc^). For 
very long cylinders of length L we can replace the sum 
over kz by an integral. The result is 



El2{<j) = 




From the argument theorem it follows that 

where /(A) is an analytic function in the complex A 
plane within the closed contour C, with simple zeros at 
Ai,A2,... within C. We use this result to replace the 
sum over n in Eq.Q by a contour integral 

X ^ In 0(A). (5) 

Here the function Q(A) is the ratio Q(A) = F{X)/F{\) 
such that F{X) vanishes at A„ and F{X) vanishes at A„ 
for all n. 

To proceed we must choose an adequate contour for the 
integration in the complex plane (see for details). We 
find 

(6) 
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As we will see, Q{iy) is a real function hence, the integral 
over y in Eq. ^ is restricted to y > k^- After some 
straightforward steps one can re-write this equation as 



El2 — 



L 
in 



dy y\nQ{iy) 



(7) 



As we have already mentioned, this expression is valid for 
conductors of arbitrary shape, as long as there is transla- 
tional invariance along the z-axis. The role of the point- 
matching method will be to provide an explicit expression 
for the function Q. 

A general solution of the scalar Helmholtz equation ^ 
can be written as: 



R{r,e) 



oo 

E 

m— — oo 



[AnJ,n(Ar)+S,„ff^i)(Ar)]e^' 



(8) 



(1) 



where (r, 6) are the polar coordinates, and Jm and 
are the m-th order Bessel functions. The constants Am 
and Bm are determined by the boundary conditions. For 
TM modes, the function R must verify Dirichlet bound- 
ary conditions on each conductor Ci and C2 (for TE 
modes, one should impose Neumann boundary condi- 
tions). The key point is to impose the boundary condi- 
tions on a finite number of points, as it is esquematically 
shown in FigO Therefore, for TM modes we have 



m— — 5 
S 

0= [^m^m(Ar,) + S„HW(Ar,)]e™''', (10) 

m=-S 

where (rp, Op) are points of the curve Ci and (r^, 6q) be- 
long to the curve C2 . We assume that with a finite num- 
ber of terms, say 25-1-1, we can acquire a desired compu- 
tational accuracy in the solution. Thus, we should satify 
Eq.® and Eq.dTU]) &t 2S + I points on Ci and C2. 



FIG. 2: Sections of the waveguides, indicating the points 
where we impose boundary conditions on each conductor 

After imposing boundary conditions, we can write, in 
matrix form, the set of linear equations for the coefficients 
Am and S,„ in Eqs. dH) and ^ 



MiA + M2B 
NiA + N2B 



where {Mi) 
H^m\\rp)e'"''- , iNi),„ 



pin 



0, 
0, 

J™(Ar,)e""«' 



(11) 



and 



{N2), 



For this system to have 



non trivial solutions, we shall ask that the determinant 
be zero. Thus, assuming that all matrices are square (i.e 
there are 25' + 1 points on each curve). 



det 



Ah 



M2 
N2 



detN2 ■ det{Mi - AhN:^^ Ni) = . 

(12) 

This equation determines the eigenfrequencies associated 
to the geometry, and the usual approach in classical elec- 
tromagnetism would be to find numerically its roots. 
However, as already mentioned, in order to compute the 
Casimir energy it is not necessary to find each eigenvalue. 
It is far more efficient to use the argument theorem, and 
obtain the sum over all eigenvalues as an integral in the 
complex plane. As this is in general a divergent quantity, 
we will compute the interaction zero-point energy by sub- 
stracting the energy of the same geometric configuration 
but with a very large outer conductor. This means that 
we will have another set of equations similar to Eas. fTT]) 
but evaluated on a different outer surface (in which 
Tq is replaced by ar^, with a ^ 1) 

det ^Jl « det ^ det[M2Nr] = . 

(13) 

In the equation above we have taken into account that, 
when using the argument theorem, all matrices will be 
evaluated on the imaginary axis. In this case, the differ- 
ent matrices become proportional to the modified Bessel 
functions 

2 

Jniiiyrp) = i"'Irniyrp); H^\iyrp) = -{~i)"'Km{yrp) . 

TT 

(14) 

After these substitutions, it is easy to check that the ma- 
trix vanishes in the limit a — > 00. Then, in order to 
compute the finite interaction energy with the argument 
theorem, the relevant determinant is 



det[{Mi - M2N^^Ni)M2\N^)-^] 



(15) 



Neglecting the self-energies of each configuration (that do 
not contribute to the interaction energy between the con- 
ductors), we obtain the expression for the function Q to 
be included in the expression for the Casimir interaction 
energy in Eq. ([7]): 

Q{iy) = det[I~AhN^^N2M2^] 

= det[I - MiN^^N2M2^], (16) 



where 



(Mi)p„ = /„(2/rp)e™«- 

(M2)p„ = Km{yrp)e'"''- 

{Ni)qm = /rn(2/r,)e""«' 

iN2)qm = A™(yr,)e™«' 



(17) 



Note that the factors that multiply the modified Bessel 
functions in Eq. fTi)) cancel out in the function Q. 
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There are similar expressions for the TE modes, in 
which each Bessel function is replaced by its derivative. 
The full Casimir interaction energy is given by the sum 
of the TE and TM contributions. 

As a first example, we consider the case of having a 
configuration of concentric cylinders. In that case, the 
points on curve Ci arc all on a circle of radii Vp ~ a, i.e. 
we arc choosing points of coordinate (a, Op). Similarly, 
points on the outer cylinder of radii h will be (6, 6q) since 
we can choose points with the same polar angle Op but 
of different radial distance. In such a case, the matrices 
can be factorized as 

{Mijpm = Qpmlmiya), {M2)pm = QpmKmiya), 
(^l)pm = QpnJmiyb), {N2)p,n = QpmKm{yb) .(18) 

Here Q is the matrix that contains the angular contribu- 
tion 



Qpm = exp{im6'p}, 



(19) 



where, in the case of circular sections, we have 6p = 
2t:p/{2S+ 1). Inserting Eqs.dH]) and (HH) into Eq.^ 
we find 



Q{iy) = det 



1 



I{vci)K{yb) 
I{yb)K{ya) 



(20) 



where we introduced diagonal matrices / and K. In 
this way, we reobtain the known expression for the TM 
Casimir energy for concentric cylinders Q : 



-.TM 



12 



L 
in 



dy y \niY\_ 



1 



lm{ya)Km{yb) 

{yb)Km(ya) 



For the TE modes we the Casimir energy is given by a 
similar expression, in which each Bessel function is re- 
placed by its derivative. 

We have developed a numerical Fortran routine in or- 
der to evaluate the Casimir interaction energy for a two- 
conductor waveguide of arbitrary cross sections. With 
a two dimensional mesh and the use of 25* -I- 1 points 
on each curve, we define the cross section of the waveg- 
uide. Once the points are selected, we impose Dirichlet 
and Neumann boundary conditions on each of them and 
define the M's and TV's matrices. Then we use a stan- 
dard Fortran routine to invert the matrices and define 
the matrix whose determinant gives the Q function (see 
Eg. pH)) ) to be included in Eq.d?]). Finally, we calculate 
the determinant and perform an integration over all val- 
ues of y. The parameters used by the program are: the 
number of points used to define the geometry 25" -I- 1, the 
integration limit (j/max) and the precision desired. 

In principle, the energy can be computed accurately if 
enough points are taken. This quantity rules the bigger 
order of the Bessel functions involved in our simulations, 
and this means that, sometimes, special precautions have 
to be taken to maintain accuracy of computation (such 
as the approximation for small argument for the modi- 
fied Bessel functions). However, we should remark that. 



though its limitation in the number of points selected, 
this technique is of great use for evaluating the TE and 
TM contributions to the energy. In the case of smooth 
geometries, we shall see that the number of points se- 
lected is enough to reproduce results in excellent agree- 
ment with previous calculations. 

Finally, wc would like to mention that the choice of 
the functions used to described the general solution to 
the Helmholtz equation (Bessel and complex exponential 
functions in our case) depends on the cross section of the 
waveguide. Other geometries could also be worked out 
using a different set of complete functions. 



III. TESTING THE METHOD: ECCENTRIC 
CYLINDERS 

We shall begin by applying the point-matching ap- 
proach to the case of two eccentric cylinders. In this 
case, we can check our simulations against our previous 
numerical work pl| . Therein, we numerically evaluated 
the Casimir interaction energy for two eccentric cylinders 
using the formula 



Ei2 = ^ I dy y 



An 



ln(M™(y))+ln(Af^™(2/)) 



TM/ 



where M™{y) 
det[S„p '^^i 



det[(5„ 



I] and M™(y) 



ATM] 

■"np ' ^n,p\ 

np — ^^.p]. The matrices A^-J^ and A™ were ana- 
lytically derived in a previous work p^ . Herein, we shall 
evaluate the Casimir interaction energy by the use of the 
point-matching technique. We will follow the notation of 
Ref.[l3l, denoting by a and b the radii of the inner and 
outer cylinders, respectively, and by e the eccentricity of 
the configuration (distance between the axes of the cylin- 
ders). The adimensional quantitites 5 = e/a and a = b/a 
will be useful to describe the numerical results. 

The comparison between the new and old approaches 
for the evaluation of the Casimir interaction energy be- 
tween eccentric cylinders was done for several runs with 
the same parameters in both programs. In this case we 
used a mesh of 21 points for simulating each cylinder, 
and matrices of 21x21 in our old approach. In this way, 
we are able to reobtain Figs. 2 and 3 of our previous 
work 2l| as can be seen in Figs. [3] and [HI where we plot 
the interaction energy for eccentric cylinders for differ- 
ent values of the radii and the eccentricity, using both 
methods. It is easy to note in Figs. [Hini and [5] that the 
difference among them is less than the 2% in the case 
of a > 2 and 5 < 0.5. The error can be blamed to the 
bigger quantity of algebraic operations contained in the 
new approach. It is worth mentioning that adding one 
more point to the point-matching method, i.e. having 23 
points in each circle, represents a variation in our result 
of only 0.000004% for a = 3 and S = 0.1 and 0.0004% 
for a = 2 and S = 0.1. Thus, we see the great agree- 
ment between both approaches. As can be seen from 



5 



3 - i \ - 3 - ■ 

i - i 6=0.1 ^ .- . 6=0.5 
2 ^y'-l ;■- 2 ]■■;:'■[ r^-,--; 
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-3 > -3 - i > \ > \ - 

-3-2-10123 -3-2-10123 

FIG. 3: Numerical mesh and points used in the point match- 
ing method for simulating the boundaries between two eccen- 
tric cylinders. In this case, we show two different values of a 
and of the eccentricity 5. For simplicity, for each point on the 
inner cylinder we choose a point on the outer cylinder with 
the same angular coordinate. 
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FIG. 5: Comparison between the new and old approaches for 
the evaluation of the Casimir energy differences as a function 
of Of = b/a for different values of the eccentricity S. Energies 
are measured in units of L/Atto^. 
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FIG. 4: Comparison between the point-matching method (set 
in the plot as new) and the numerical method used in Ref. [21l | 
(set as old) for numerically evaluating the Casimir energy as 
a function of the eccentricity 5 for different values of q = 6/a. 
Ai5i2 refers to the energy difference between the eccentric and 
the concentric configurations. Energies are measured in units 
of L/47ra^ 



Figsin]and[71 small deviations appear for large eccentric- 
ities and small a. We expect the precision in these cases 
to be improved by considering grids with a larger num- 
ber and/or a non- uniform distribution of points on the 
surfaces. 

In Fig. [S] we show the comparison between the new 
and old approaches for concentric {6 = 0) cylinders, as a 
function of a. In this case, the agreement between the 
full numerical and the point-matching approaches is even 
better. 
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FIG. 6: Ratio between the Casimir energies evaluated using 
the new and the old approaches, as a function of the eccen- 
tricity S for different values of a — b/a. 



IV. CYLINDRICAL RACK AND PINION 

When two concentric cylinders have corrugations, the 
vacuum energy produces a torque that could, in princi- 
ple, make one cylinder rotate with respect to the other. 
This "cy lindrical rack and pinion" has been proposed in 
Ref. [28| , where the torque has been computed using the 
proximity force approximation. It was further analyzed 
in [2^, where the authors obtained pcrturbative results 
for Dirichlet boundary conditions in the limit of small 
amplitude corrugations. In this Section, we numerically 
evaluate the Casimir interaction energy for two concen- 
tric corrugated cylinders. The cylinders have radii a and 
&, and we will denote by r_ — b — a the mean distance 
between them and by r_|_ ~ a + b the sum of the radii. As 
in the previous Section, we will use the notation a = b/a. 
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0.98 
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FIG. 7: Ratio between the Casimir energies evaluated using 
the new and the old approaches, as a function of a = 6/a for 
three different values of S. Energies are measured in units of 
L/Aiva'. 



which reads 



El2 

nr+L 



cos(z^0o) 



240r° 



h^Bl'^^'iy), (22) 



where y ~ r-/r^ and B^^^ {y) is given by 

BP'^iy) = ^ E / ^dx 

^yi — — 



m— — oo 
2 



Ay 



1 



1 



(1 - 2/2) Dm{y,x) D„i+„iy,x)' 

The functions Dm are given by 

Dm{y;x) = /„,(a;[l +y])i^,„(.x[l -7/]) 
- -y])i^,„(.T[l +7/]). 



(23) 



(24) 
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2 4 6 8 10 12 14 16 18 20 
a 

FIG. 8: Comparison between the new and old approach for 
numerically evaluating the Casimir interaction energy in the 
concentric case {S — 0), as a function of a = b/a. The point- 
matching approach is greatly accurate in this case. 

In Fig. Owe show two geometries with different frequen- 
cies associated to the corrugations: on the left side, v = i 
and on the right side = 5, both for a = 2. The points 
in the mesh are described by the following functions: 

ha{e) hsui{ve) 

hb{9) = hsm{i^e + (f>o), (21) 

where h is the corrugation amplitude and ly is the fre- 
quency associated with these corrugations. The Casimir 
torque can be calculated by taking the derivative of 
the interaction energy with respect to the shifted angle 
T = —dEi2/d4>Q. In the case of the TM mode (Dirichlet 
boundary conditions), details of the perturbative calcu- 
lation can be found in Ref.jl^. Therein, the authors 
obtained an analytical expression for the Casimir inter- 
action energy as a function of the angle (jy^) for small /i, 




-3-2-10123 -3-2-10123 



FIG. 9: Numerical mesh and points used in the point match- 
ing method for simulating the boundaries between two con- 
centric cylinders with harmonic corrugations. In this case, we 
show a = 2, 4>o = 0, for two different values of v. 

In Figliniwe show the numerical evaluation of the TM 
Casimir interaction energy for this geometry. The plot 
shows the results obtained using our point matching ap- 
proach with a = 2 and corrugation frequency v = for 
different values of the amplitude of the corrugation h. 
As expected the amplitude of the oscillations grows with 
h. For each value of h we have performed a numerical 
fit of the data in order to compare with the analytical 
prediction presented in Eg. ip^ . With dotted lines we 
have plotted the fit y{x) = A * cos{x) for each curve in 
FigHni The agreement between dots and dotted lines is 
extremely good. 

In Table U we can sec the comparison between the an- 
alytical prediction and our numerical results for different 
values of a, u and h = h/a. For each simulation we have 
used a mesh of 37-41 points, depending the value of h. 
For a = 2 and h = 0.1, having 31 or 37 points selected 
to define each boundary has a relative error of 0.00004% 
in the final value of the energy. For i/ = 5, this difference 
is approximately 0.002%, as one can anticipated (for a 
smoother curve less points are needed in order to achieve 
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FIG. 10: Casimir interaction energy (TM modes) as a func- 
tion of 00 for a = 2 and different values of perturbation 
h = h/a. The different shaped dots are the numerical data 
obtained with our programe while the line represents the nu- 
merical fit of each curve. Energies are measured in units of 
L/a^. 



V 


a 


h 


A (Analytical) 


A (Numerical) 


3 


2 


0.01 


0.000030414 


0.00003044 


3 


2 


0.05 


0.000760 


0.00077 


3 


2 


0.1 


0.0030 


0.0033 


3 


2 


0.3 


0.027 


0.078 


3 


3.5 


0.01 


0.000003925 


0.000003928 


3 


3.5 


0.05 


0.0000981 


0.0000993 


3 


3.5 


0.1 


0.000392 


0.000412 


3 


3.5 


0.3 


0.00353 


0.00567 


5 


2 


0.01 


0.00002046 


0.00002049 


5 


2 


0.05 


0.000511 


0.000528 


5 


2 


0.1 


0.00204 


0.00232 


5 


2 


0.3 


0.018 


0.071 



TABLE I: Comparison between the analytical fEq. (|22|l ') and 
numerical predictions (fit of the form y(x) = A* cos(a::)) of 
the Casimir interaction energy (Dirichlet modes) for different 
configurations of the concentric corrugated cylinders. 



It is worth emphasizing that there are no analytical pre- 
dictions for this mode so far, being this the first evalua- 
tion of the TE Casimir interaction energy for corrugated 
concentric cyhnders. 
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FIG. 11: Casimir interaction energy (TE modes) as a function 
of 00 for a value a = 2 and different values of perturbation 
h. The different shaped dots are the numerical data obtained 
with our programe while the line represents the numerical fit 
of each curve. Energies are measured in units of L/a^ . 
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the same accuracy in the result). In the table, we can 
see that for small values of h, the values of the amplitude 
of the numerical fit of the data are extremely similar to 
those predicted analytically. However, we must note that 
this technique does not have a constraint on the value of 
h so far. For bigger values of h one has to include more 
points in the mesh so as to have a well defined geometry 
and mantain accuracy in the results, that differ from the 
analyical predictions. 



In addition, in Fig. [TT] we have presented the evalua- 
tion of the Casimir interaction energy for the Neumann 
(TE) modes. Therein, we see that the behaviour is qual- 
itatively similar to that of the Dirichlet modes, with a 
different value for each numerical fit of the data points. 



FIG. 12: Casimir interaction energy (TE and TM modes) as 
a function of 0o for a = 2, = 3 and h = 0.3. The dif- 
ferent shaped dots are the numerical data obtained by our 
programme while the line represents the numerical fit of each 
curve. In this case, the plot shows that the exact result can- 
not be fitted by a function y{x) = A * cos(2;). Energies are 
measured in units of L/a? . 

It is worth to remark that, when the amplitude of the 
corrugation is not very small, the exact results cannot 
be reproduced with a simple fit of the form y{x) = A ^ 
cos(a:). This is illustrated in FigfT2l where we see that, 
for the biggest corrugated amplitude that we included in 
Table|l](/i = 0.3), the exact result differs from the cosine 
function. 
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FIG. 13: TM and TE contributions to the Casimir interaction 
energy in the concentric corrugated case. Dirichlet (TM) con- 
tribution is bigger than the Neumann (TE) one (in absolute 
value) . 



Finally, in Fig[T3] we show the different contributions 
of the TM and TE modes to the interaction energy for 
bigger values of a. As expected from previous results 
for large values of a, Dirichlet contribution is bigger 
(in absolute value) than the Neumann one. 



V. OUTER CONDUCTORS WITH FOCAL 
LINES 

Some time ago, there was a conjeture [3^ based on 
a geometric optics approximation, about the possibility 
of focusing vacuum fluctuations in parabolic mirrors. It 
was argued that a parabolic mirror is capable of focusing 
the vacuum modes of the quantized electromagnetic field, 
therefore creating large physical effects near the mirror's 
focus. The physical manifestation of this focusing is a 
growth in the energy density and mean-squared electric 
field as the focus is approached. In particular, the en- 
ergy density would diverge as the inverse fourth power of 
the distance from the focus in the case of perfect conduc- 
tivity. These results would imply that the focused vac- 
uum fluctuations will enhance Casimir forces on atoms 
or other particles near the focus. The sign of the force 
could draw particles into the vecinity of the focus [lOl . 

With this motivation, in this Section we shall evalu- 
ate the Casimir interaction energy for configurations in 
which the outer conducting shell has a cross section that 
contains focal points. 



denote by a the radius of the inner cylinder, by bi and 62 
the minor and major semiaxes of the ellipse, respectively, 
and by / the distance between the foci and the center of 
the ellipse. The coordinates of the center of the cylinder 
with respect to the center of the ellipse will be (ex,ey)- 
We will use an additional tilde to denote adimensional 
quantities, i.e distances in units of a: hi ~ bi/a , f = f /a, 
etc. 

For this configuration, we use a mesh like the one pre- 
sented in Fig. I14i where we show an inner cylinder, and 
an outer ellipse with semiaxes bi ~ A and 62 = 4.33. The 
ellipse has two focal points at / = 1.66. We present the 
results for the Casimir energy in Figs ll5l and 1161 
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FIG. 14: Mesh used to evaluate the boundaries condition in 
the framework of the point-matching method. We represent 
a small centered cylinder and vary the position of the outer 
ellipse. Parameters used: minor ellipse semiaxis 61 = 4, major 
ellipse semiaxis 62 = 4.33 and focal position / = 1.66. 
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A. Cylinder inside an ellipse 

To begin with, we will compute the Casimir interaction 
energy between one small inner cylinder and an outer 
ellipse, by the use of the point matching method. We will 



FIG. 15: Numerical evaluation of the Casimir interaction en- 
ergy for an inner cylinder an eccentric outer ellipse, as a func- 
tion of the position of the cylinder along the vertical axis. 
Energies are measured in units of L/a^ . 
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FIG. 16: Moving the inner cylinder along the minor semiaxis 
the ellipse we show that the center of the ellipse is an unstable 
equilibrium position. This is shown for different positions of 
the inner cylinder along the major semiaxis. Energies are 
normalized to the value of \Ei2{tx = 0)|. 

From Figlin] it is possible to see that there is an un- 
stable equilibrium position at the origin under displace- 
ments of the inner cylinder along the (vertical) ey direc- 
tion. As expected, it is also possible to check that the 
energy grows as well as the cylinder gets closer to the 
surface of the outer ellipse. Fig[T5] also shows a mono- 
tonic behaviour of the energy as a function of the posi- 
tion, even when passing through the focus. Otherwise, 
Fig ll6l shows the unstable equilibrium position at the ori- 
gin when moving the inner cylinder in the (horizontal) 
direction. It is also important to stress that, when con- 
sidering horizontal displacements at a fixed vertical po- 
sition, the higher the altitude of the inner cylinder, the 
narrow the inverted potential in Fig ll6l 

B. Cylinder inside a Parabola 

Following the same idea than in the previous sub- 
section, here we describe the numerical computation of 
the Casimir interaction energy for a small inner cylinder 
with circular cross section, inside a large cylinder with 
parabolic section. In Fig llTl wc show the two dimensional 
cross section of the mesh used to simulate this geometry 
and impose the boundary conditions. The parameters for 
this case are the radius of the inner cylinder, a, and the 
focal distance of the parabola, /. As before we introduce 
f = f/a. 

In Fig llSi we show the behaviour of the Casimir in- 
teraction energy with the position of the inner cylinder. 
The number of points used is approximately 40. It is 
easy to see that the bigger increasement of the energy 
(and therefore of the force) is in the direction labeled as 
ex- We can also note an increasement in the energy as the 
inner cylinder gets closer to the vertex of the parabola (in 
the direction), a biproduct of the proximity between 



FIG. 17: Mesh used to evaluate the boundaries condition in 
the framework of the point-matching method. We represent 
a small centered cylinder and vary the value of the position 
of the outer parabola. Parameter used / — 4. 

the inner cylinder surface and the one corresponding to 
the parabola. On the other hand, as the behaviour of 
the energy is symmetric in the ey direction, the force Fy 
vanishes on the horizontal axis. There is also an unstable 
equilibrium position at a particular point on this axis, as 
it is suggested by FiglTHl As in the previous example, 
the energy and the force do not have a special behaviour 
near the focus. 

-0.0005 I 1 , 1 , , , , 1 , 1 

-0.001 - - 

-0.0015 - 



-0.002 - 

-0.0025 - 

-0.003 - 

-0.0035 - 




-0.004 - , \ i 

-0.0045 - / I 

-0.005 I ' ' ' ' ' ' ' ^ ' 

-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

e 

FIG. 18: Numerical evaluation of the Casimir interaction en- 
ergy between an inner cylinder inside a bigger parabola. We 
move the inner cylinder both in directions and ey . Param- 
eter used: f — 4. Energies are measured in units of L/a^. 



VI. FINAL REMARKS 

In this paper we presented a new numerical method to 
compute the vacuum energy for arbitrary geometries with 
translational invariance. The method is based on the use 
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of the point-matching approach, in which the boundary 
conditions are imposed on a discrete set of points on the 
surfaces of the conducting bodies. This approach is com- 
bined with the argument theorem, in order to trade the 
sum over eigenvalues for an integral in the complex plane. 

After testing our method against previous results, we 
have computed the Casimir interaction energy for new 
geometries. In the case of the cylindrical rack and pinion 
described in Section IV, we have seen that, for corruga- 
tions of small amplitude, the numerical results for the en- 
ergy show a harmonic dependence with the shifted angle, 
and a quadratic dependence with the amplitude of corru- 
gations, in agreement with previous analytic perturbative 
evaluations. This behaviour disappears for larger ampli- 
tudes, where the exact results show more pronounced 
peaks. 

In Section V we computed the Casimir interaction en- 
ergy between an inner cylinder and outer surfaces with 
focal lines. The motivation for looking at these config- 
urations was to see whether there was a non trivial be- 
haviour of the energy near the focus or not. For both 
cases considered (ellipses and parabolas), we have found 
that the energy and forces are monotonic across the focal 
lines. This may be a peculiarity of the geometries with 
translational invariance considered here. We have also 
confirmed the existence of unstable equilibrium positions 
of the inner cylinder, that coincide with the location sug- 
gested by simple geometric arguments. 

The examples discussed in this paper illustrate the 
simplicity and power of this approach. We have used 
a straightforward version of the point-matching method, 
with a naive choice of the points on the curves (note that 
in all examples we have chosen pair of points with the 
same angular coordinate with respect to the inner cylin- 
der). For less symmetric configurations, and when the 
surfaces of both conductors are closer to each other, it 
will be necessary to consider grids with a larger number 
of points, and to optimize their positions. As in the appli- 



cations to acoustic or classical electromagnetism, special 
care must be taken for surfaces with pronounced edges, 
clefts or " handles" , where the point-matching technique 
may not be accurate to determine the eigenfrequencies. 

It is certainly possible to go beyond the geometries 
considered here. Still considering geometries with trans- 
lation invariance along the z-axis, it would be possible to 
analyze waveguides with more than two conducting bod- 
ies. The generalization for non-perfect mirrors is also 
possible. For the case of scalar fields satisfying arbitrary 
matching conditions on the surfaces (rather than Dirich- 
let or Neumann boundary conditions), the generalization 
is relatively straightforward. However, for the electro- 
magnetic field the calculation is more cumbersome, since 
in general it will be not possible to treat independently 
the TE and TM modes. 

The point-matching approach can also be generalized 
to three dimensional compact objects. The starting point 
in this case would be the solutions of the three dimen- 
sional Helmholtz equation, written in terms of an ade- 
quate basis (products of spherical Bessel functions and 
spherical harmonics). We expect these geometries to re- 
quire much more computational effort, and a more so- 
phisticated method to optimize the choice of the grid of 
points on which the boundary conditions are imposed. 
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